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g ABSTRACT 

0957+561 is the first gravitational lens system to be discovered and 
first for which a time delay was measured. Because the system is unusu- 
ally rich in observables such as image positions, fluxes, VLBI structure, 
and polarization it has been modelled quite extensively. However, since 
it resides in a cluster providing substantial convergence, the velocity dis- 
persion of the lensing galaxy is an important component to models. We 
have measured this dispersion to be 288 ±9 km s -1 (1-sigma), which im- 
plies a Hubble constant of 72 ±7 km s _1 Mpc -1 (1-sigma) according to the 
preferred FGS model of Grogin and Narayan, and 70 ± 7kms" 1 Mpc" 1 
(1-sigma) using the SPLS model. The unknown velocity dispersion 
anisotropy of the central galaxy produces an additional uncertainty of 
perhaps 15% or more. We see no variation in dispersion with position 
to a radius of ±3". In addition to a galaxy from the background cluster 
at z = 0.5, our slit serendiptously picked up a galaxy at z = 0.448, and 
two galaxies which may be part of a background cluster at a redshift of 
z = 0.91. 



1 Based on observations at the W. M. Keck Observatory, which is operated jointly by the California Institute of 
Technology and the University of California 
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Subject headings: cosmology — distance scale — gravitational lensing - 
quasars: individual (0957+561) 
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1. Introduction 

Gravitational lenses bend light from a back- 
ground source and can cause multiple images 
of the source to appear. Light has a differ- 
ent time of flight along each path, and if the 
source should vary there will be a time de- 
lay between when each of the images changes 
brightness. Such a time delay can transform 
dimensionless redshifts and angles into phys- 
ical distances, hence provide a Hubble con- 
stant (Refsdal 1964). In principle this is an 
extremely powerful and important method for 
measuring the cosmological parameters em- 
bodied in distances and redshifts, party be- 
cause of the potential accuracy at high red- 
shift but mainly because they suffer a com- 
pletely different set of biases from more tra- 
ditional measures of H . 

A time delay for the lensed system 0957+561 
has been known for a decade (Schild and 
Cholfin 1986), although Press, Rybicki, and 
Hewitt (1992) disputed their value and only 
recent measurements by Kundic et al. (1997) 
have settled the issue beyond doubt. See 
Haarsma et al. (1997) for a review of the 
time delay measurement. Models of the sys- 
tem have also undergone improvement; Gro- 
gin & Narayan (1996) [GN] have carried out a 
particularly extensive analysis which has been 
widely used for estimating H . Romanowsky 
& Kochanek (1998) have recently examined 
the question of how velocity anisotropy in the 
lensing galaxy affects the use of a central ve- 
locity dispersion to estimate the lensing po- 
tential. 

The models of GN can yield a Hubble con- 
stant either by using an estimate of the lo- 
cal convergence k due to the cluster poten- 
tial or by measuring the depth of the lens- 
ing galaxy potential from its velocity disper- 



sion a. The first measurement of the veloc- 
ity dispersion was by Rhee (1990) who found 
300 ±50 km s^ 1 . Recently, observations with 
LRIS on the Keck telescope by Falco et al. 
(1997, FSMD) have found a dispersion which 
drops rapidly with radius, from 316 ± 14 at 
0.2" from the center of the lensing galaxy Gl, 
dropping to 266 ± 12 at 0.4". Using the lat- 
ter value they derive a Hubble constant of 
Hq = 62 ± 7kms~ 1 Mpc~ 1 , and they spec- 
ulated that a massive central object might be 
causing an anomalously large central disper- 
sion. 

The very rapid drop in velocity dispersion 
with distance is unusual behavior for such a 
massive galaxy, however, which inspired an 
attempt to reobserve the lensing galaxy in 
0957. One concern about the observations 
of FSMD centers on their choice to run the 
slit through the Gl galaxy and the QSO B 
component simultaneously. While this offers 
both the galaxy and QSO light in one expo- 
sure it suffers from much more scattered light 
and could potentially corrupt the galaxy spec- 
trum in unforseen ways. Indeed, FWMD sug- 
gest redoing their observation with the spec- 
trograph slit rotated by 90 degrees. 

We chose to do just that: rotate the slit 
perpendicular to the line connecting Gl and 
the QSO B component, excluding as much 
QSO light as possible. In addition, we at- 
tempted to match the QSO contribution leak- 
ing into the slit as accurately as possible by 
taking alternate exposures with the slit ro- 
tated by 180° and offset to the other side of 
the QSO B component, thus obtaining a rela- 
tively pure QSO spectrum which illuminates 
the slit in the same way as during the Gl ex- 
posure. 
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2. Observations and Reductions 

0957+561 was observed on March 30 and 
March 31, 1997 using the Low Resolution 
Imaging Spectrograph (LRIS) (Oke et al. 1995) 
at the Keck II telescope on Mauna Kea, along 
with calibration exposures, an observation of 
the cluster MS1358+62 to act as a veloc- 
ity dispersion calibrator, and HD 132737 and 
AGK2+ 14873 as radial velocity templates. 
The observations are summarized in Table 
1. The sky was clear and the seeing was 
about 0.8" throughout both nights. Long 
slits of 1.0" and 0.7" were used along with 
the 600 1/mm grating blazed at 5000A. The 
600 1/mm grating was rotated either to cover 
4800-7680A for the QSO and cluster obser- 
vations or 3800-6080A for the template stars. 
The spectral resolution was 4.65A fwhm for 
the 600 1/mm grating with a 1.0" slit, 3.59A 
fwhm for the 600 1/mm grating with a 0.7" 
slit, and the scale along the slit was 0.211"/pixel. 
The template stars were guided smoothly across 
the slit in many locations in order to have uni- 
form illumination across the slit, to build up 
signal to noise, and to map out loci of con- 
stant slit position across the detector. 

Figure 1, derived from a very deep image 
by Fischer et al. (1997), illustrates where the 
three slit positions lay. PA108 runs through 
the lensing galaxy, Gl, two neighboring galax- 
ies, t2 and t3, and also catches two other 
galaxies at the end of the slit, tl and t4. Since 
we needed as clean a spectrum of Gl as possi- 
ble we rotated the slit by 180° to PA288 and 
offset to a point symmetrically located on the 
other side of QSO image B. This also yielded 
a spectrum of another galaxy, t5 as well as tl 
again. Finally, we offset to center the slit on 
QSO image A to acquire a pure QSO image, 
and we also picked up galaxies t6, t7, and t8. 




Fig. 1. — Illustration of the slit positions and 
identifications of the galaxies observed here. 
North is up and East is left; G2-G5 are the 
designations of Young et al. (1980); and t6 is 
galaxy 20 and t2 is galaxy 21 from Angonin- 
Willaime et al. (1994). 

The spectra were reduced using software 
described in detail by Tonry (1984). The ba- 
sic steps are to flatten the images, remove 
cosmic rays, derive a wavelength solution as 
a function of both row and column using sky 
lines (wavelengths tabulated by Osterbrock et 
al. 1996), derive a slit position solution as a 
function of both row and column using the po- 
sitions of the template star images in the slit, 
rebin the entire image to coordinates of log 
wavelength and slit position, add images, and 
then sky subtract. A quadratic fit to patches 
of sky on either side of the object (including a 
patch between for the MS1358 observations) 
did a very good job of removing the sky lines 
from the spectra. 

As has been stressed by Franx (1993) and 
by van Dokkum and Franx (1996), measur- 
ing the velocity dispersion of a galaxy at a 
redshift significantly larger than zero must 
be done carefully, since the instrumental res- 
olution of a spectrograph tends to a con- 
stant number of angstroms, whereas the red- 
shifted spectrum has been stretched. Hence 
a straightforward cross-correlation or fitting 
method will underestimate the dispersion of 
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the galaxy. Kelson et al. (1997) measured dis- 
persions in the cluster MS 1358+62 using high 
resolution template spectra and very careful 
modelling of the spectrograph resolution de- 
rived by measuring sky line widths. The ob- 
servations here use the trick that the ratio 
between an 0.7" and 1.0" slit is slightly larger 
than (1 + z) for these lenses, hence a template 
measured with the 0.7" slit will have almost 
exactly the same instrumental resolution as a 
galaxy at a redshift of z ps 0.3. This is borne 
out by the ratio of the measured spectral res- 
olutions: 4.65A + 3.59A = 1.30. 

As a test that this observing procedure will 
give correct dispersions, we also observed the 
cluster MS1358+62 for which Kelson et al. 
(1997) measure a dispersion of 303 ±6 for the 
central galaxy, Gl, and 219 ± 4 for the neigh- 
boring bright galaxy, G2. We discuss our ob- 
servations of MS1358 in Tonry (1998), but 
the salient point is that the velocity disper- 
sions measured for the MS 1358 galaxies using 
the stellar templates reveals no significant off- 
set from Kelson's values, Gl: 299 ± 22; G2: 
235 ± 23. We therefore conclude that use of 
the stellar templates is likely to give us accu- 
rate velocity dispersions. 

2.1. Companion Galaxies 

Extraction of the spectra of the companion 
galaxies required no correction for QSO light. 
In each case the spectrum was extracted, and 
cross-correlated with the template spectrum 
according to Tonry and Davis (1979) as well 
as being analyzed by the Fourier quotient 
method of Sargent et al. (1977). For each 
spectrum the redshift, error, and velocity dis- 
persion were calculated. The redshift calcu- 
lation included the entire spectrum, whereas 
the spectrum blueward of 4000A (rest frame) 
was excised for the dispersion calculation, 



since the calcium H and K lines, intrinsically 
broad and subject to interstellar absorption, 
are always problematic for dispersions. Dis- 
persions are not corrected for any aperture 
effects, hence correspond to an aperture of 
approximately 1" square. 

Three galaxies, t4, t7, and t8, showed un- 
mistakable emission lines at 7126, 7176, and 
7105A respectively, and t7 also had a strong 
emission line at 5396A, as well as ones which 
might be present at 7037A and 7248A, al- 
though they are confused with night sky lines. 
We believe that the emission in t4 is [Oil] 
3727A b ecause of the corroboration from ab- 
sorption lines at the red end of the spectrum 
which correspond well to the UV lines be- 
tween [Oil] and Ca K. Galaxies t7 and t8 
only received an integration of 500 sec, so it 
is more difficult to be certain of their redshift. 
The four emission lines in t7 correspond well 
to [Oil] at 5376A and Hp and the [OIII] lines 
in the red, so we are quite confident about 
the redshift we assign it. Galaxy t8 is much 
more problematic, with only one emission line 
present. We tentatively identify this line with 
[Oil], giving it a redshift similar to that of 
galaxy t4, but a deeper integration will be 
necessary before we have much confidence in 
this redshift. 

Table 2 lists the redshifts, errors, velocity 
dispersions, errors, and cross-correlation sig- 
nificance "r" values for each spectrum. The 
spectra are illustrated in Figure 2. Compar- 
ison of the redshifts here with those mea- 
sured by Garrett et al. (1992) and Angonin- 
Willaime et al. (1994) shows complete agree- 
ment within the errors. 

2.2. Lensing Galaxy 

Anticipating a serious contamination of the 
spectrum of the lensing galaxy Gl by QSO 
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Fig. 2. — The spectra are shown for the four 
galaxies at redshift 0.355, t6, t5, t3, and t2 
(offset by 400, 250, 150, and for clarity), 
and the four galaxies at higher redshift, tl, 
t4, t7, and t8 (offset by 800, 400, 150, and 0). 
These spectra have been Gaussian smoothed 
with a fwhm of 6A. The top axis shows rest 
wavelength at a redshift of 0.355. 
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Fig. 3. — The combined spectrum of the lens- 
ing galaxy Gl and the leakage of light from 
QSO component B are shown, as well as the 
scaled spectrum of QSO component A and the 
difference which is the "pure" spectrum of Gl 
which is subjected to redshift and dispersion 
analysis. These spectra have been Gaussian 
smoothed with a fwhm of 6A. The top axis 
shows rest wavelength at a redshift of 0.355. 

component B, we had collected light from the 
other side of B, with the slit arranged so that 
the light gradient across the slit would match 
the exposures of Gl. However, the seeing was 
good enough that there was surprisingly lit- 
tle contamination (less than 30 percent, and 
the QSO light is mostly featureless), and on 
the night of Mar 31 we dispensed with the 
matching exposures of B. A simple scaling of 
the exposure of QSO component A removed 
the narrow absorption lines of Mgll and gave 
us the "pure" spectrum of Gl illustrated in 
Figure 3. The spectrum has been averaged 
over 11 rows, equivalent to 2.3 arcsec. The 
effective aperture is therefore 1 arcsec by 2.3 
arcsec. 

We used three different methods for com- 
puting the velocity dispersion of the lensing 
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galaxy Gl. The first is the standard analy- 
sis which uses the template star observations. 
The second method uses the MS 1358 galaxies 
as velocity dispersion templates. The third 
method involves fitting the sky lines to deter- 
mine the instrumental profile and then creat- 
ing a synthetic template from higher resolu- 
tion data. 

Using a star as a velocity dispersion tem- 
plate requires a narrower slit to make the in- 
strumental resultion commensurate for tem- 
plate and galaxy. As demonstrated by the 
comparison between the dispersions measured 
for MS1358-G1 and MS1358-G2 with the higher 
resolution data of Kelson et al., our trick of 
using an 0.7" and 1.0" slit is likely to be pretty 
good at the redshift of 0957, but we expect 
that there may be systematic differences in 
dispersion according to which wavelengths are 
compared. We found no difference between 
the four observations of Gl, so report only 
the average dispersion which we find to be 
290 ± 20. Using the blue half of the template 
gave a dispersion of 317 and the red half gave 
285, so we view this result with some caution. 

The second method exploits the MS1358 
data again, cross-correlating the lensing galaxy 
against each of MS1358-G1 and MS1358-G2. 
Then, cross-correlation of MS1358-G1 against 
MS1358-G2 with varying amounts of extra 
broadening allows us to calibrate the width of 
the correlation peak of each as a function of 
velocity dispersion since we know the disper- 
sions of each of these galaxies. The advantage 
here is that 0957 is at a small redshift relative 
to MS1358, so the difference in instrumental 
resolution is negligible. The disadvantage is 
that this result uses low signal-to-noise tem- 
plates and depends on the accuracy of Kel- 
son's dispersions. Using MS1358-G1 as a tem- 
plate, we find a dispersion of 285 ±20 for 0957- 
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Wavelength (A) 

Fig. 4. — The normalized spectrum, the fit 
derived from the Fourier fit (offset by —0.15), 
and the difference (offset by +0.5) are plotted 
as a function of observed wavelength. The big 
feature at 5840A is the G band. The quality 
of the fit is clearly very good. 

Gl, and using MS1358-G2 as the template, we 
find a dispersion of 283 ± 23. 

Finally, we use the method introduced by 
van Dokkum and Franx (1996) and measure 
the width of sky lines, using them to broaden 
a synthetic, redshifted, high resolution stellar 
template which we provide to the Fourier fit- 
ting method of Franx et al. (1989). Several 
template stars were fitted to the galaxy spec- 
trum, and only the best fitting template was 
used to derive the final velocity dispersion, 
and we fitted over the wavelength region of 
5558- 6710A to avoid the atmospheric B band 
at 6800A. This reduction also reproduces Kel- 
son's dispersions of MS 1358 very accurately 
and gives a velocity dispersion of 289 ± 9 for 
0957-G1. The spectrum and the fit are shown 
in Figure 4. The quality of the data and the 
fit are clearly excellent. We regard this as the 
most reliable method, but since the other two 
methods are consistent, we offer a weighted 
average of 288 ±9 as our best estimate for the 
velocity dispersion of the gravitational lens. 
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Fig. 5. — The velocity dispersion is plotted 
as a function of position along the slit (open 
symbols are results from use of stellar tem- 
plates; star symbols show results from line fit- 
ting). The data points from FSMD are also 
shown as small, filled symbols. 

The dispersion as a function of position 
along the slit is illustrated in Figure 5. Us- 
ing the stellar templates, we find values of 
286±30, 340 + 20, 294 + 10, 293±10, 287+10, 
308 ± 15, and 291 ± 30 at positions of -3.2", 
-1.8", -0.7", 0.0", +0.7", +1.8", and +3.2" 
along the slit, and using the sky line calibra- 
tion we find values of 311 ± 28, 293 ± 14, 
286 + 11, 301 ± 14, and 248 ± 24 at posi- 
tions of -1.7", -0.6", 0.0", +0.6", and +1.7", 
i.e. no significant variation of dispersion with 
radius. The discrepancy between the two re- 
duction methods is reasonably consistent with 
the formal errors, and the two points at ±3" 
should not be taken very seriously. 

3. Discussion 

Our redshift for galaxy tl agrees well with 
the measurement by Garrett et al. (1992), 
confirming the existence of the cluster they 
found at z — 0.5. If the redshift of galaxy 
t8 is confirmed to be 0.907, it suggests that 
galaxies t4 and t8 might be part of a back- 



ground cluster along the line of sight. The 
spectrograph slit was about 160" long and we 
would probably detect a galaxy within ±l"of 
the slit, so the total area surveyed by the three 
slits was about 500 square arcsec. The angu- 
lar size distance at z — 0.91 is 1355 Mpc 
(for Q = 0.25), hence 2 galaxies at that 
distance projected within 0.02 Mpc 2 , is a re- 
markably high density. The rest frame veloc- 
ity dispersion suggested by the two galaxies 
is about 800 km s -1 . 

The velocity dispersion derived here differs 
somewhat from that reported by FSMD. As 
shown in Figure 5, they found a dispersion 
which dropped very rapidly from 316 km s^ 1 
at 0.2" = OMh^kpc (for fi = 0.25) to 
266 km s" 1 at 0.4" = 1.7 kpc. This would 
be extremely unusual behavior for a large el- 
liptical galaxy, not to mention a spectrum 
which has been convolved with the seeing and 
instrumental profile. We do not confirm their 
rapid drop in a, and do not agree with their 
conclusion that the lensing galaxy harbors a 
massive, dark object. We believe therefore 
that their results are subject to substantial 
systematic error (perhaps because of scat- 
tered light), that the error that they assign 
to their dispersion measurement is much too 
optimistic, and that our measurement is prob- 
ably more reliable. 

Grogin and Narayan (1996) offer two for- 
mulae for calculation of H from a time de- 
lay and a velocity dispersion. Using the 
time delay of 1.142 ± 0.004 year (1-sigma) 
from Kundic et al. (1997), and a veloc- 
ity dispersion of 288 ± 9, the two models 
give 72 ± 7kms" 1 Mpc- 1 (FGS model) and 
70 ± 7 km s- 1 Mpc" 1 (SPLS model). The 
velocity dispersion and model give roughly 
equal contributions to the error. It would 
be quite possible to improve the accuracy 
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of the velocity dispersion by reobserving the 
lens at higher dispersion, perhaps with a 
1200 line/mm grating, and it is also possi- 
ble that the models need improvement. For 
example, Bernstein et al. (1997) use HST 
WFPC2 images of 0957 to measure the center 
of the lensing galaxy, and find that its loca- 
tion is inconsistent with the GN model. 

The main uncertainty in the modeling is 
due to anisotropics in the velocity disper- 
sion. Richstone and Tremaine (1984) showed 
that the measurement of velocity dispersions 
through a central aperture does not constrain 
the mass of an elliptical galaxy very well. 
From their figures, we estimate the uncer- 
tainty in the mass at 40% for our aperture 
of 1 x 2.3 arcsec and an effective radius of 4.5 
arcsec (Bernstein, Tyson, & Kochanek, 1993). 

A recent analysis by Romanowsky & Kochanek 
(1998) has confirmed this additional uncer- 
tainty. They find that the mass cannot be 
well measured if only a central velocity dis- 
persion is measured. They show that an ac- 
curate measurement of the lineprofiles can 
constrain the mass much better, but unfor- 
tunately, no such measurements are avail- 
able for Q0957. Our data give uncertainties 
on the order of 0.04 on the Gauss-Hermite 
terms used for such an analysis (e.g., van der 
Marel & Franx 1993). When Romanowsky 
& Kochanek impose the constraint that the 
Gauss-Hermite terms do not deviate from 
mean Gauss-Hermite terms measured for lo- 
cal galaxies, they find that the uncertainty in 
H is 011 the order of 15 %. Given these uncer- 
tainties, it is clear that single lenses will not 
give a very accurate value for H , but larger 
samples may provide a good mean. 

Despite these uncertainties, it is hearten- 
ing to find that we are beginning to converge 
on a value for the Hubble constant via very 



different paths. Recent values which depend 
on the Cepheid calibration include 82 ± 8 
(Lauer et al. 1997, SBF of HST images of 
BCGs), 81 ± 6 (Tonry et al. 1997, SBF of 
nearby galaxies), 73 ± 6 ± 8 (Freedman 1997, 
reporting on the Cepheid calibration of a va- 
riety of secondary distance estimators includ- 
ing supernovae and Tully-Fisher) , and 69 ± 5 
(Giovanelli et al. 1997, Tully-Fisher relation). 
Other values which are nominally indepen- 
dent of the Cepheid calibration include 64±20 
(Schechter et al. 1997, PG1115 time delay), 
and 54 ± 14 (Myers et al. 1997, Sunyaev- 
Zeldovich effect). This is not the forum to 
discuss the relative merits of various methods, 
but it is clear measurements of the Hubble 
constant are finally attaining some degree of 
maturity, and that measurements from gravi- 
tational lens time delays can be an important 
contributor to the quest for H . 

We are grateful to P. Fischer for provid- 
ing us with the exceedingly deep image of the 
0957 field used in Figure 1. As always, we are 
deeply appreciative of the capabilities of the 
Keck telescopes and the LRIS spectrograph. 
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Table 1 
Observing Log. 



Objects 


UT Date 


UT 


sec 


z 


PA 


Exposure 


Slit 


AGK2+ 14783 


3/30 


5:39 


1 


,01 


90 




1.0 


AGK2+ 14783 


3/30 


5:46 


1 


.01 


90 




0.7 


0957 GL 


3/30 


6:14 


1 


,31 


108 


1500 


1.0 


0957 QSO B 


3/30 


6:54 


1 


,26 


288 


1500 


1.0 


0957 GL 


3/30 


7:28 


1 


.24 


108 


1500 


1.0 


0957 QSO B 


3/30 


7:59 


1 


.24 


288 


1500 


1.0 


0957 QSO A 


3/30 


8:26 


1 


.24 


288 


500 


1.0 


MS1358 G1+... 


3/30 


11:12 


1 


.37 


24 


1500 


1.0 


MS1358 G1+... 


3/30 


11:39 


1 


.36 


24 


1500 


1.0 


HD132737 


3/30 


15:23 






90 




0.7 


HD132737 


3/30 


15:27 






90 




1.0 


AGK2+ 14783 


3/31 


5:14 


1 


.00 


90 




1.0 


AGK2+ 14783 


3/31 


5:26 


1 


.01 


90 




0.7 


0957 GL 


3/31 


6:55 


1 


.25 


108 


1500 


1.0 


0957 GL 


3/31 


7:25 


1 


.24 


108 


1500 


1.0 



Note. — PA is east from north, exposures are in seconds, slit widths 
are in arcseconds. 
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Table 2 
Redshifts and Dispersions. 



ID 


Y# 


G# 


A# 


-<) 




z 


± 


a 


± 


r 


emission 


Gl 


88 






-0.3 


1.1 


0.3562 


0.0001 


294 


10 


12.6 




tl 


40 


R33 




60.9 


20.4 


0.5044 


0.0005 


193 


70 


4.7 




t2 


100 




21 


-11.2 


-2.6 


0.3550 


0.0001 


105 


66 


14.6 




t3 


115 






-22.5 


-6.3 


0.3583 


0.0001 


72 


100 


14.6 




t4 


164 






-79.0 


-24.7 


0.9120 


0.0003 








[Oil] 


t5 


116 


R19 




-22.8 


-8.3 


0.3532 


0.0001 


154 


63 


8.3 




t6 


97 


R13 


20 


-7.5 


3.1 


0.3547 


0.0002 


68 


143 


6.9 




t7 








21.6 


8.9 


0.4476 


0.0004 








[OII],^,[OIII] 


t8 








27.4 


13.9 


0.9070? 


0.0005 








[Oil]? 



Note. — Columns: Galaxy name, Identification numbers from Young et al. (1981), Garrett 
et al. (1992), and Angonin-Willaime et al. (1994), Position with respect to QSO component B 
(arcsec), redshift and error, velocity dispersion (kms -1 ) and error, cross-correlation r value, and 
emission lines observed. 
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